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Abstract 

We provide elements for a discussion of the expected v signal from Supernova Rem- 
nants (SNR) in the Milky Way. After recalling why SNR are interesting and cer- 
tain remarkable achievements of H.E.S.S., we describe a simple and straightforward 
method to evaluate the v fluxes from 7-ray data. For an ideal detector, we get a flux 
of 5 thoroughgoing muons per km 2 per year from RX J1713. 7-3946 in ANTARES 
location and above E t h = 50 GeV; similar calculations for Vela Jr show that the 
number of events, to be evaluated precisely after the next detailed observations of 
H.E.S.S., are larger. We comment on the role of neutrino oscillations. 

1 The cosmic ray/SNR connection 

Supernovae are suspected to be the cosmic ray (CR) accelerators since '34 
(Baade & Zwicky pQ). 30 years later, Ginzburg & Syrovatsky 2 a remarked 
that if 10 % or so of the SNR kinetic energy £$n ~ 10 51 erg (=1 foe, also 
known as 1 bethe) goes in CR, the Milky Way losses are compensated: 

VCR PCR q 1 x ^ SN 

where Vcr = itR 2 H (R = 15 kpc, H = 5 kpc), tqr = W 7 yr and t$n = 
30 yr. The 'diffusive shock wave acceleration' mechanism based on Fermi 
ideas jH] is being developed to explain CR acceleration in SNR (lecture of 



1 



Blasi); CR accelerate in an expanding shock wave of size R = u t (it ~ 

5.000 km/s), mostly active in the first f,000 yr, determined by M e j ec t a ~ 

R 3 n ISM - 

Hillas selected a list of open (connected?) questions [3]: 

• How to "inject" electrons? ( diffusive shock accel. is incomplete? ) 

• Why isotropy? How T = 2.1 — > 2.7? ( propagation/reacceleration? ) 

• .Emax? ( R ~ Dsohm/u limit, countered by Bell & Lucek |3J/ ) 

• Too few point sources of VHE 7? ( something lacking? } 

• How to firmly exclude a leptonic origin? ( TeV v are needed? ) 

We will be mostly concerned with the last questions (the last one is not in 
Hillas's list, but perhaps he would have included it, if giving this talk). 

1.1 The landscape after the first H.E.S.S. results 

For a few young shell-type SNR observed in VHE 7, the "hadronic" hy- 
pothesis seems plausible; more crucial tests will be possible with future 
observations by H.E.S.S. (VERITAS, MAGIC) and other instruments. The 
closest SNR's, whose properties we recall here, are of particular interest: 
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1° 
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no(t yet?) 
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2 kpc 


36' 


1,000 yr 


Cas A 


HEGRA (maybe) 


58°08' 


3 kpc 


6' 


320 yr 



(note however that the "distance" and "age" are not reliably determined, 
see e.g. [§]). For the first two the velocity of expansion is as expected, for 
the other it is a bit larger. RX J1713. 7-39346 and Vela Jr are the most 
intense shell-type SNR in the TeV 7 sky; thus, it is particularly important 
to discuss the expected neutrino fluxes from them. 

The best known 7 spectrum is the one of RX J1713. 7-3946. H.E.S.S. [7j 
showed that the spectrum of RX J1713. 7-3946 deviates from a power law 
distribution above ~ 10 TeV. Assuming that the 7-rays come from CR 
interactions with the environment (are of hadronic origin), the cut in the 
CR (proton) spectrum should be around 150 TeV, in agreement with naive 
expectations from diffusive shock acceleration models. Specific models of 
this SNR have been proposed: Malkov, Diamond, Sagdeev '05 [S] suggest 
that the nearby molecular cloud has a main role for CR interactions whereas 
Berezhko & Volk '06 9 fit H.E.S.S. observations starting from the opposite 
view; in both models the hadronic contribution is sizable or dominant. In 
short, H.E.S.S. observations renewed the interest in the CR/SNR connection 
and perhaps, figured out the most intense sources of VHE 7 rays. 



Two web links that can help to keep information updated are: 

* H.E.S.S. Source Catal., www.mpi-hd.mpg.de/hfm/HESS/ [W. Hofmann]; 

★ Catalogue of SNR, www.mrao.cam.ac.uk/surveys/snrs/ [D. Green]. 
See also the review of H.J. Volk ^U] and the lecture of Wei Cui. 



Motivated by the (shell-type, young) SNR / CR connection and by the exist- 
ing plans for large neutrino telescopes, we calculate the flux of TeV neutrinos 
from the SNR with known VHE j-ray spectrum. Indeed, during CR accel- 
eration the SNR are transparent to their 7 radiation. Thus, we can convert 
the measured 7 ray flux (from ir° and rj) into an expectation for the neu- 
trino flux (from 7r ± and K^) under the hypothesis that the radiation is of 
hadronic origin. We begin by discussing flavor oscillations, describe the 7/1/ 
connection and estimate the rate of events in (km 2 class, ideal) neutrino 
telescopes. 

2.1 Oscillations 

The flux of neutrinos-from meson decays-are modified by the oscillations: 



The oscillation probabilities take the simplest form, Gribov-Pontecorvo's 
(namely, the one that applies for low energy solar neutrinos): 



There is no MSW effect ^2j, for matter term is negligible close to the SNR 
and too large in the Earth. With central values of the mixing elements 17% 
we get Pun ~ 0.4 and P efl ~ 0.2; that is, 1 /2 of the original z/„ and Vn fluxes 
reach the detector. 

We performed also a detailed (or sophisticated) analysis 



2 TeV neutrinos from SNR 
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Pw = E l^l with £ > * = e >v,T 



i=l 



£(P wt ) oc max e 



(iW- p w( fl )) 2 



x C osc .(6) with a — ► 



where 9 are the measured parameters taken from |13j that we marginalize 
away by maximizing the result. We get P^ = 0.39 ± 0.05 and P efl = 
0.22 =p 0.05 where most of the error (0.04) is due to #23. 



To understand the uncertainty budget one can use an expansion in the 
small parameters |14j : 



Pnn - V 2 - V 8 - y and - V 4 + y, 

-- sin 2 20i2, 

= COS 2# 2 3 x /4 + 013 COS 5qP V x ( 1-a: )/2. 



where 



In view of astrophysical uncertainties and small counting rates we believe 
that the uncertainties in the oscillation parameters do not have an important 
role for the discussion of SNR v and presumably even for other cosmic 



sources. 1 



2.2 The connection between 7 and v 

For RX J1713. 7-3946 there are various calculations in the literature: 

1. Alvarez-Muhiz h Halzen '02 [IB] inspired by CANGAROO first obser- 
vations use F 7 oc E~ 2 and obtain F v = F® oc F 7 by 



dE u E v F v (Ey) = J dE'y E~/Fj(E~i 

;mm/12 JE™ in /6 



2. Costantini & V '04 Jl] use F-y oc E 2 ' 2 as extrapolated from early 
H.E.S.S. below 10 TeV results and adopt standard techniques (see 

e.g., Ha to 



AX 2Z r , 
A,, 



-F 7 = — ^—F p and similarly for F v 



Both methods, however, are tailored for power law spectra; and we know 
from the only detailed observation that we have (RX J1713. 7-3946) that this 
is not a good approximation. Thus we are lead to recalculate the neutrino 

x We comment on a proposal ^5] to study 613 and Sep with a source of v £ ^jj H7| 
using the v flux ratio F M / (F e + F T ) = P eAl /(l — Pep). Since the shift due to 5cp and 86iz 
around 0° is SP efi — 0.02, first we should know the impact of 623. If this were negligible, 
the number of e + r signal events N should obey SN/N < 10 %. This needs 60 years of 
an ideal detector and a small systematic 8b /b < 1 %, if we use the event rates per year 
in a km 2 area of |15| s = 16 (signal) and b = 145 (background): 10 3 signal events over 
10 4 events. We get instead s — s e + s T — 1.3 (s e = N t j^ a rTeV a{E) F e (E) dE; targets: 
Nt = 4.5 10 38 nucleons; flux from US], fig.l: 7 10 6 (E/TeV)" 3 1 1/TeV yr km 2 with 
P ee w 0.6; cross section: 4 10 _36 (i?/TeV) ' 87 cm 2 ). Even assuming the source is real, 
with s — 1.3 the question is whether a signal can be seen. 



fluxes. In principle, one could de-convolute the CR flux from the 7-ray 
flux and then obtain the neutrino flux, as described in the lecture of K.- 
H. Kampert. In our case, when the atmosphere (the target for CR) is 
much thinner-it is transparent-a much simpler and direct approach |2J is 
possible based on the classical techniques of Lipari '88 |22| . In fact, from 
the integral expression for VHE 7-rays F 1 (E) = J™ dE' 2 F 7T o(E')/E' we 
find immediately: 

Due to the approximate isospin-invariant distribution of pions, F n = F^o ~ 
F w + ~ E w -, we find for the neutrinos from n + — > y^v^: 

JE/{X-r) l-r E' 2(1 - r) 

where r = (m^/m^) 2 . The neutrinos v = v„, u e from /x + decay are: 

F„(E V ) = [ ^ F^iE,,) [g (y) - P^) gi(y)} where ^ = — (3) 
Jo y y 

Qi are polynomials, F^ and P^ (polarization averaged over tt distribution) 
are also known. The contributions of semileptonic decays to the z/s 
and the one of r\ — > 77 to 7-rays can be described with a contribution 
proportional to the pionic one; these two correct the formulae in opposite 
directions so that effectively one should just add a contribution of kaons of 
the order of 10 percent. For details, see [2"T] . 



2.3 A case study: RX J1713.7-3946 

The neutrino flux, evaluated as described above, is shown in figure ^ and 
tabled For the purpose of illustration, we use it to calculate the number of 
muon events Nu, + Nn for an ideal detector 2 using: 

POO 

N„ = A-T-f liv - / dE u F Vfl (E u )Y^E u , E th )(l - a v .{E u )) 

where E v is the neutrino energy before the interaction point and the various 
quantities in the previous formula are defined as follows: 

2 By definition, 'ideal detector' means that all muons above threshold are detected: 
e(_E M ) = 1 for Efx > E t h- Of course in a real detector (where the geometrical area A 
changes with the time of observation t and, above all, the efficiency e increases with the 
energy) the impact of the cut in the spectrum is expected to be much more important. 



Figure 1: spectra, corresponding to two fits of the H.E.S.S. VHE j-rays from RX 
J 17 13. 7-3946 J^: a broken power law and a power law with exponential cutoff; the 
third fit- a curved power law-is incompatible with a hadronic origin, for it increases 
before 40 GeV. The flux is the same in our approximation. 
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• A = 1 km 2 and T = 1 solar year. 

• The source is below ANTARES horizon (=visible) for fu v = 78 %. 

• The threshold for muon detection is E t h = 50 GeV (as low as possible). 

• The muon range (that goes in the yield Y^) is calculated for water. 

• The neutrino absorption coefficient a v , averaged over the daily location 
of the source, is calculated for standard rock. 

We find that the total number of events does not depend crucially on 
the extrapolation: 

jy _ f 4.8 per km 2 per year [exponential cutoff] 
M ^ | 5.4 per km 2 per year [broken power law] 

This can be compared with the 9 events in ^1] (power law extending till 1 
PeV) and the 40 events in [18. (oscillations, livetime and absorption ignored). 
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Table 1: Differential flux from RX 31713.7-39^6 for selected energies (1 st and 
4 lines). Units: TeV for the energy, le-12/TeV s cm 2 for the flux. 



3 Summary and perspectives 



We discussed the expected neutrino flux from 7-transparent accelerators of 
cosmic rays, and emphasized the case of young SNR. We showed that for 
RX J1713.7-3946 (the best known SNR in TeV sky, thanks to H.E.S.S.) the 
expectations are stable: ~ 5 events per km 2 per year in an ideal detector. 
The median neutrino energy is 3 TeV. Since the detected jx are softer than 
the (i in the production point (that in turn is softer than the impinging v) 
several events will fall in an energy region where the atmospheric background 
and the role of imperfect detection efficiency are important: see the lectures 
of Lipari and Lucarelli. Thus, we believe that it would be desirable to have 
a detailed discussion of the characteristics of a detector that aims to see a 
v signal from RX J1713.7-3946. 

Sometimes soon H.E.S.S. should tell us more on another intense VHE 
7-ray source, RX J0852. 0-4652 (Vela Jr). This SNR has F 7 = 6.5 million 
7-rays i£ _2,1 /(TeV km 2 yr) below 10 TeV. If the exponential cutoff is at 
-E^cut = 50 (150) TeV (and again, if these 7 rays are of hadronic origin) we 
get Nn + Nn = 10 (14)/(km 2 yr) in an ideal detector, with a significantly 
higher energy |2Jl . 
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A Details of signal evaluation 

Following ^1] and we describe here two specific elements used for the 
evaluation of the signal and comment on them. 

A.l Muon background 

It should be possible to increase the acceptance by designing an angular cut 
that depends on the energy and on the time of the event. Yet, here we follow 
the simplest prescription to exclude cosmic muon contamination: accept 
only the events below the horizon of the detector. The angle 9 between an 



astronomical object and the vertical of the detector is 

cos9(t) = ( cos <5, 0, sin <5 ) • ( cos 0cos(7rt/r), cos <ft sm(irt / t) , sin (ft ) 

where 1) 5 is the declination of the object (for the two sources discussed 
in the text: RX J1713.7-3946 = -39°46', RX J0852.0-4652 / Vela Jr 
= — 46°52'); 2) (ft is the latitude of the detector (for the detectors in 
the Northern hemisphere: Baikal = 51°50', ANTARES = 42°50', NEMO 
= 36°30', NESTOR = 37°33'); 3) 2 x r = 23^56 m 4 s is the duration of the 
sidereal day; 4) t is the time measured from the point when the object is at 
the apex. If cos#(i = 0) < 0, the object remains always below the horizon; 
if the converse happens, but cos 8{t = r) < 0, the object is observable for a 
fraction of the time 1 — Tq/t, determined by the condition cos # (to) = 0. 

A. 2 Earth absorption 

The Earth absorption of the neutrino depends on the time t through the 
average column density in the Earth and it is mostly due to CC interac- 
tions. NC interactions are ~ 1/3 smaller and, furthermore, cannot remove 
a neutrino, but only degrade its energy. A detailed evaluation of this effect 
is numerically demanding [2Sj, however it supports the expectations from 
a simple heuristic argument inspired to the formalism of 'scaling' (proper 
of strong interaction): NC interactions increase the absorption, but only by 
few %. Thus in a first approximation NC effect can be neglected. 

A. 3 Comments 

In the light of the new H.E.S.S. results, it is useful to note that: 

1. Accepting events from a few degree above the horizon (that is, setting 
the cut to a certain value cos 9q > 0) can lead to significant increase of 
the events, especially for RX J1713. 7-3946 [23j. E.g., for ANTARES 
the fraction of time becomes fu v = 88 % just adding 5 degrees above 
the horizon. 3 

2. The absorption coefficient obtained by the approximate procedure 
mentioned in sect. I A. 2l using a power law distribution E~ 2 ' 2 and aver- 
aged over the fraction of the time when cos8(t) < is very similar for 
RX J1713. 7-3946 and for Vela Jr (slightly weaker in the second case). 

3 For a cosmic muon that reaches the depth of 1 km under the surface, 5° means more 
than 12 km of water but recall that a muon can change its direction by scattering. 



3. If the neutrinos are confined to be of relatively low energy, as it is the 
case of RX J1713. 7-3946, the effect of Earth absorption is small. 
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Discussion 

AURELIO GRILLO: 1) How reliable is the 7-transparency hypothesis? 
2) What is the relative weight of oscillations and correct spectrum? 

FRANCESCO VISSANI: 1) For a SNR as RX J1713.7-3946 it is pos- 
sible to check that the matter is very diffuse either distributing several solar 
masses in the wide volume, or more directly using X-ray data, so that the 
absorption of 7-rays is negligible. E.g., for the molecular cloud: 100 Mq at 
1 kpc and a size of 15' means a column density of about 2 10 20 protons/cm 2 
(10 protons per cm 3 ). For other sources as /iQSO (namely, stellar black 
holes or neutron stars with jets) this does not apply and a relatively weak 
7-radiation could be compatible with intense v fluxes |25j . 4 
2) Again for RX J1713. 7-3946, each of the two effects reduce the signal in an 
ideal detector by a factor of 1/2. However oscillations are universal whereas 
deviations from power-law depend on the individual object. We hope this 
is less pronounced for Vela Jr, and wait for the response of H.E.S.S. 

GIANNI NAVARRA: Is it possible to "extrapolate" these calculations 
in order to obtain a neutrino luminosity from the galactic disk, in the hy- 
pothesis that all CR are produced by SNR? 

FRANCESCO VISSANI: Strictly speaking the results in eas ing are 
just a tool to assess an upper bound to v from 7-transparent sources (or if 
you like it more, a lower bound to v from VHE 7 ray sources of hadronic 
origin) and I can only hope that they will be useful for tasks as the one you 

4 It was suggested that fj, (but not 7T ) are absorbed in /iQSO. This leads to a peculiar 
flavor neutrino ratio, as for v from hypothetical neutron sources (see footnote^ . Of course, 
before testing flavor ratios we should be sure that we can at least observe muon neutrinos 
above TeV and the most cautious attitude would suggest to perform all possible tests with 
conventional means first. 



propose. 5 One ingredient for the calculation seems to be the following: if we 
have a new SN each 30 years and if they are mostly active for some thousand 
years, we have 40-80 SNR that inject CR effectively in each moment. These 
are most likely located in the spiral arms of the Milky Way, but we can get a 
rough idea of their distribution assuming that their density is r exp(— r/ro), 
with ro ~ 3 kpc |2S]- In this way we can find the expected distance of the 
closest one, their average distance, etc. These results could/should be cross- 
checked with SNR databases. An obstruction could arise in the description 
of an 'average' SNR: at different ages the intensity and the distribution of 
7-rays could change. In fact, I suspect that the SNR that we are beginning 
in the 7-rays sky are peculiar objects, as those associated with molecular 
clouds, and/or the brightest /closest ones. 

BERND ASCHENBACH: Just a comment: Even if it should turn out 
that the TeV spectrum in RX J1713. 7-3946 is dominated by hadronic pro- 
cesses, it is also very useful to compute the contribution of leptonic processes, 
since it should be there if adiabatic diffusive acceleration (involving protons 
and electrons) at the SNR shock front is universal. The question is at which 
level. 

FRANCESCO VISSANI: Certainly this task is of paramount impor- 
tance. Despite the difficulties, it should profit from the accumulation of 
new observations, it can be approached in concrete models of CR acceler- 
ation (such as [H| and [§]) and it is presented in a moderately optimistic 
perspective in a recent authoritative review work |32j . 

But I would like to use this occasion to recall that, beside the contribu- 
tion from leptons to 7 rays spectrum, there is also a hadronic contribution to 
the spectrum of hard X-rays. This originates from the non-thermal popula- 
tion of electrons produced in muon decays (incidentally, no new formula 
is needed to describe the flux of electrons or positrons for this is just the 
same as the one of v e or v e , see eq. EJ) that could provide us with a further 
handle to disentangle the components of leptonic and hadronic origin. 



5 Just after Vulcano 2006 several papers appeared addressing this type of problem 
|27l 1281 129 j . mostly using assumptions on the CR flux at least for intermediate steps. 
More information and discussion in two subsequent conferences 1301 and |31| . 



